ACS Catalysis

In Anabaena variabilis (Trichormus variabilis) phenylalanine ammonia-lyase (AvPAL), a conserved lid-like loop sits over the active site and has been studied both for its role in positioning a catalytic tyrosine and for its contribution to phenylalanine aminomutase (PAM) activity. While the active site architecture and substrate specificity of AvPAL have been extensively characterized, the dynamic behavior of this unstructured loop beyond its role in catalysis remains poorly understood. Here, we investigate the functional role of this loop by restricting its mobility through targeted interchain disulfide bond engineering. Three in-house approaches were designed to predict ideal cysteine residue pairs: (i) quantifying pair interaction energies via electrostatic and van der Waals forces, (ii) generating a contact map of residues within 5 [A] proximity, and (iii) implementing a machine-learning model trained on datasets from PDBCYS, SPX, and an internal database to rank cysteine pair likelihood within disulfide bond geometric constraints. Our machine-learning-guided strategy yielded a successful variant with complete oxidation efficiency in E. coli. Rigidification of this loop reveals that it also functions as a regulator of substrate specificity. Multi-scale molecular simulation analyses (molecular dynamics, metadynamics, quantum/molecular mechanics) reveal that this modification alters the active-site pocket by reducing the conformational dynamics of substrate binding. Our findings underscore the delicate balance between enzyme flexibility and catalytic efficiency, providing novel insights into the role of this understudied dynamic loop region in AvPAL.

OXA-232, an OXA-48 like carbapenemase stands amongst newly identified beta-lactamases that causes of the extensive of beta-lactam resistance. While active-site residues are well characterised, the contributions of conserved non-active-site residues in exerting enzymatic activity remain unexplored, limiting our understanding about the roles of these residues in the overall OXA-232 function. To address these gaps, the conserved residues S118, V120, L158, and D159 of OXA-232 positioned adjacent to the active-site motifs and within the omega-like loop were substituted with alanine. Substitutions of S118A and D159A rendered the expressing cells susceptible to penicillins, cephalosporins, and carbapenems, whereas the cells harbouring OXA-232V120A and OXA-232L158A proteins exhibited substrate-selective susceptibility changes. Kinetic analysis with purified proteins revealed the reduction in catalytic efficiency of all the mutants compared to wild-type protein. Though the L158A and D159A mutated proteins become deacylation-deficient, the mutations S118A and V120A exhibited selective acylation defects without trapping intermediates. It is evident from circular dichroism spectroscopy and molecular dynamics simulations that OXA-232S118A, OXA-232V120A, and OXA-232L158A nearly retained their secondary structures and compactness, except for OXA-232D159A, which presumably triggered a misfolding leading to destabilisation of the omega-loop. Interestingly, bicarbonate supplementation partially rescued the lost activities in soluble mutants, underscoring the carbamylation dependence. Taken together, these findings establish S118 and D159 as essential for core catalysis and structural integrity, with V120 and L158 modulating substrate-specific turnover and orientation. The current study reappraised the mechanistic insights of OXA-48-like carbapenemases, providing significant resources in rationally designing future therapeutics to combat carbapenem resistance.

Genetic code expansion introduces new-to-nature chemical moieties into ribosomally synthesized proteins. In practice, the scope of functional groups that can be accessed using this method is often limited by noncanonical amino acid (ncAA) availability. Producing ncAAs directly in cells can circumvent poor ncAA uptake or commercial unavailability, but limited enzymes suitable for this application exist. In vitro evolution campaigns have been remarkably successful in yielding synthetically useful "ncAA synthases." However, these enzymes are optimized for preparative-scale synthesis and their activities often do not translate well to cellular biosynthesis. Thus, expanding strategies to engineer enzymes specifically for ncAA production within cells will benefit further implementation of genetic code expansion. Here, we use phage-assisted noncontinuous and continuous evolution to evolve enzymes for improved synthesis of non-canonical tyrosine derivatives in E. coli. Using simple serial passaging, we uncovered mutations that doubled the production of an expensive ncAA, 3-methoxytyrosine, by tyrosine phenol lyase, and furthermore evolved variants that enable 3-iodotyrosine biosynthesis, a transformation the parent enzyme is unable to catalyze. Additionally, we evolved a recently reported tyrosine synthase for improved production of 3-halogenated tyrosines, identifying variants that exhibit high activity even at low substrate concentrations owing to a [~]8-fold reduction in KM. Our results demonstrate that phage assisted evolution can be used to rapidly improve the activity of enzymes for ncAA production in cells.

Thiol dioxygenases (TDOs) catalyze the incorporation of molecular oxygen into thiol metabolites and N-terminal cysteine residues of regulatory proteins, thereby playing critical roles in sulfur metabolism and oxygen sensing. Despite extensive study over the past two decades, the molecular basis for substrate recognition and the catalytic mechanism of TDOs remains controversial, owing to the scarcity of substrate-bound structures and direct evidence for catalytic intermediates. Herein, we present a comprehensive study of mercaptosuccinate dioxygenase (MSDO), a TDO originally identified in Variovorax paradoxus B4, using a combination of structural, biochemical, spectroscopic, and computational approaches. MSDO oxidizes both (S)- and (R)-mercaptosuccinate (MS) with similar Km values but exhibits approximately 2.5-fold higher turnover for the (S)-enantiomer. Crystal structures of MSDO reveal that both (S)- and (R)-MS coordinate the iron in a bidentate mode via their thiolate and proximal carboxylate groups, with the distal carboxylate adopting distinct orientations. Two active-site Arg residues recognize the substrate carboxylate groups and thereby stabilize a flexible C-terminal loop, underpinning a catalytic site gating mechanism in MSDO. EPR spectroscopy corroborates bidentate coordination, showing conversion of a high-spin {FeNO}7 complex to a low-spin species upon substrate binding. Time-resolved in crystallo reactions capture two key iron-bound intermediates, namely an unprecedented monooxygenated sulfenate and a dioxygenated sulfinate product. These structural snapshots are supported by DFT calculations that point to a stepwise oxygen atom transfer pathway. Computational analysis further accounts for the kinetic differences between the substrate enantiomers, as rationalized by structural comparisons, active-site geometry, and second coordination sphere interactions. Together, these results elucidate fundamental principles of TDO catalysis and advance our understanding of nonheme iron-dependent oxygen activation.

O2-tolerance is a desirable property for [FeFe] hydrogenases, which are highly efficient H2-producing catalysts. While most such enzymes are highly sensitive to aerobic environments, a small number of explored representatives exhibit exceptional stability and even H2-producing activity under oxygenic conditions. However, the genetic signatures of the O2-tolerance in this class of enzymes remain largely unknown. To address this knowledge gap, we explored a close homologue of a well-characterized O2-tolerant [FeFe] hydrogenase from Clostridium beijerinckii (CbHydA1) - a hydrogenase from Terrisporobacter glycolicus (TgHydA1). Our investigation indeed confirms that TgHydA1 can transition to the O2-stable Hinact state, a hallmark of O2 tolerance. The surprising outcome is that despite the high amino acid similarity, TgHydA1 shows a substantially higher propensity to remain in the Hinact state than CbHydA1. Using protein film electrochemical experiments, we demonstrate that the root of this behavior lies in roughly tenfold slower reactivation rates than those of CbHydA1 at any applied potential. This degree and direction of variation in reactivation kinetics have not been observed before for any other O2-tolerant [FeFe] hydrogenases or their variants to date, uncovering a yet-to-be-explored facet of reactivity alteration available to these enzymes. Overall, the results presented here highlight the importance of a holistic analysis of [FeFe] hydrogenase sequences in the context of their interaction with O2 that encompasses the protein environment and properties of the auxiliary metallocofactors.

Directed evolution has long been constrained by complex screening hardware and labor-intensive workflows. Here, we report the first genuine test-tube screening platform that uses His6-tagged peptide-functionalized magnetic beads and Fe3+-decorated E. coli cells to establish a phenotype-genotype linkage, thereby decoupling ultrahigh-throughput screening from specialized instrumentation and democratizing directed evolution. The platform demonstrated a screening throughput of > 108 events s-1 and an enrichment factor of up to 63-fold. Using galactose oxidase as a model, we identified variants with up to a 26-fold increase in catalytic efficiency. Extensions to D-amino acid oxidase and alcohol oxidase yielded variants with up to 5383-fold and 25-fold improvements over their respective wildtypes after a single round of screening. These results highlight the platforms capacity to rapidly engineer H2O2-generating oxidases and to advance AI-driven enzyme design through rapid data generation.

Cytochrome P450s catalyze a diverse array of reactions including crosslinking of aromatic side chains in the biosynthesis of ribosomally synthesized and post-translationally modified peptides (RiPPs). ApyO is a cytochrome P450 enzyme that forms a C-C bond between two tyrosines in a YLY motif in the substrate ApyA, the precursor peptide of the RiPP aminopyruvatide. We utilized cell-free translation to generate ApyA variants and probe the substrate tolerance of ApyO. Through Alphafold-based modelling and in vitro assays, we show that ApyO accepts the 10 C-terminal residues of ApyA and requires a conserved Arg/Lys in the substrate peptide. Inspired by substrate sequences found in orthologous biosynthetic gene clusters, we substituted one of the tyrosine residues with a tryptophan and observed that ApyO catalyzed the formation of an N-C bond between the indole of Trp and the C{varepsilon}2 of Tyr. ApyO unexpectedly catalyzed formation of a C-O bond between the two tyrosine residues when we substituted the leucine residue in the YLY motif with tyrosine and tryptophan. We also show that a peptide containing a biaryl linkage and the C-terminal aminopyruvate displayed sub-nanomolar inhibitory activity against selected proteases. Overall, this study demonstrates plasticity in the manner of macrocyclization catalyzed by the P450 ApyO and provides a starting point for chemoenzymatic approaches towards producing diverse macrocyclic scaffolds.

Polysaccharide utilization loci (PULs) have been a goldmine for the characterization of novel carbohydrate active enzymes (CAZymes) and the understanding of their synergistic degradation of complex polysaccharides. We collected PUL predictions containing CAZymes from glycoside hydrolase families GH29, GH50 and GH117, expected to participate in marine polysaccharide breakdown. We explored the evolutionary diversity in these families in terms of sequences and PUL composition, based on sulfatases and CAZymes. From 41 selected PULs, more than 400 putative enzymes were produced, purified and screened on a large collection of carbohydrates. We attributed a function to more than 130 enzymes from five sulfatase subfamilies, 29 known CAZymes families and discovered an activity for 4 families previously of unknown function, including an -L-galactosidase structurally and functionally characterized with mutants. Finally, our detailed analysis of the enzymatic synergies in five PULs, two targeting marine polysaccharides and three targeting eukaryotic polysaccharides, by marine and human gut organisms, highlight the efficiency of our exploratory strategy.

Modern protein design methods based on deep learning allow generation of customized protein scaffolds with diverse geometries and functionalities. Here, we capitalize on these recent advances to develop hyper-thermostable de novo CO2 reductases featuring a cobalt porphyrin IX cofactor (CoPPIX). CoPPIX containing enzymes were assembled in vivo through media supplementation with cobalt salts and assessed for photocatalytic CO2 reductase activity. We identified two cysteine-ligated designs that exhibit high activity (>1000 turnovers at rates of up to 25 min-1) while suppressing competing hydrogen evolution pathways. A 2.1 [A] crystal structure shows close agreement to the design model with the Co-Cys bond programmed as intended. This study showcases the power of computational protein design in developing artificial enzymes to activate challenging molecules such as CO2.

Enzymes catalyze reactions with remarkable specificity and can unlock recalcitrant feedstocks that are dilute, complex, and variable in their constituent molecules. While characterized enzymatic reactions cover a wide range of chemistries, there are an undetermined number of cryptic activities for every known one. These cryptic activities can be elicited through rational design, adaptive laboratory evolution, and increasingly, generative models of proteins. However, prior to tuning a catalyst one must efficiently predict viable novel reactions. In this work we leverage the growing amount of mechanistic enzyme information, specifically the Mechanism and Catalytic Site Atlas, to construct a set of reaction rules that can meet this demand. By explicitly utilizing mechanistic information, the rule sets developed here more accurately identify molecular structures required for catalysis compared to existing curated and heuristically constructed rules. The 899 Distilled rules are constructed directly from characterized mechanisms and cover 62.5% of reactions from Rhea. The Learned rule set is generated from a classifier trained on mechanistic data, allowing full coverage of Rhea and precise identification of mechanism-required atoms (ROC-AUC = 0.98). Additionally, our Learned rules exhibit a more favorable tradeoff between novelty and feasibility and provide users with fine-grained control over this tradeoff. The rules are compatible with all SMARTS-based reaction network expansion and retrosynthesis software.

Caulobacter crescentus is a gram-negative bacterium that produces the anionic sphingolipid ceramide diphosphoglycerate that can substitute for the lipopolysaccharide component of the outer membrane. ccna_01210 is a gene in the operon for ceramide diphosphoglycerate synthesis and encodes for the enzyme, CpgD. CpgD is a magnesium-dependent CTP:phosphoglycerate cytidylyltransferase that catalyzes the synthesis of CDP-glycerate, and a pyrophosphate byproduct. CpgD displays substrate specificity for the nucleotide CTP and a preference for 2-phospho-D-glycerate over other phosphoglycerate enantiomers and isomers. Here, we present five high resolution structures for CpgD in various catalytic states that rationalize the specificity and preference of CpgD for its two substrates. This includes structures of apo CpgD, a CpgD-CTP-Mg2+ ternary complex, and three CpgD-CDP-glycerate-pyrophosphate-Mg2+ product bound complexes. The structures reveal CpgD nucleotide specificity is mediated by favorable hydrogen bonding interactions with the cytosine nucleobase of CTP, while the preference for 2-phospho-D-glycerate occurs due to favorable van der Waals interactions with the 2D enantiomer and unfavorable steric clashes with the 2L enantiomer. A catalytic mechanism involving a pentacoordinate transition state is proposed based on the observed stereochemical inversion of the -phosphate in the substrate CTP in comparison to the -phosphate of the product CDP-glycerate. Overall, this provides insights into the catalytic mechanism, nucleotide specificity, and enantiomeric substrate preference of the cytidylyltransferase CpgD that participates in a unique pathway of bacterial sphingolipid synthesis.

{beta}-galactosidases (BGs) are essential enzymes widely used in the food industry, particularly in the production of lactose-free products. Among them, the BG from Aspergillus oryzae is of industrial relevance due to its activity at acidic pH and moderate thermal tolerance. However, enhancing its catalytic performance remains a key challenge. Traditional enzyme engineering methods are time-consuming and resource-intensive, limiting their scalability. Recent advances in Artificial Intelligence (AI), particularly those based on Natural Language Processing, offer a promising alternative by enabling efficient exploration of protein sequence space and prediction of beneficial mutations. In this study, we introduce an ensemble-based, zero-shot Protein Language Model pipeline that reconciles predictions from six independent models (ESM2 and the five ESM1v variants) combined with a diversity-aware candidate selection strategy. Applied to the BG from A. oryzae, this approach identified beneficial mutations leading to novel enzyme variants with up to a four-fold increase in catalytic efficiency on oNPGal, a two-fold increase on lactose, and, independently, a T338I variant with markedly enhanced thermostability ({approx}80% residual activity after 24 h at 60 {degrees}C), all without requiring supervised fine-tuning on experimental fitness data. Our results demonstrate that consensus across an ensemble of PLMs can efficiently enrich beneficial substitutions in industrially relevant enzymes and substantially reduce the number of wet-lab candidates that need to be screened. Table of Contents graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=106 SRC="FIGDIR/small/726739v1_ufig1.gif" ALT="Figure 1"> View larger version (29K): org.highwire.dtl.DTLVardef@18084f7org.highwire.dtl.DTLVardef@99a102org.highwire.dtl.DTLVardef@19a64forg.highwire.dtl.DTLVardef@1f59cff_HPS_FORMAT_FIGEXP M_FIG C_FIG

Arylmalonate decarboxylase (AMDase) stereoselectively converts disubstituted malonates to chiral carboxylic acids, but its substrate spectrum is very limited regarding the size of the smaller substituent. Inspired by the observation that (S)-selective AMDase variants also convert larger substrates, we unlocked the synthesis of the (R)-enantiomers of -aryl and -alkenyl n-butanoic and n-pentanoic acids, respectively, in exquisite enantiopurity.

Photoswitchable ligands enable photocontrol of biomolecular activity by binding to targets in an isomer-dependent, light-responsive manner. Recent developments in ionizable photoswitchable ligands greatly expand their applications but introduce a major design challenge: light-responsive binding can depend on isomeric form, chemical substitution, and binding-induced shifts in protonation equilibria. These effects are tightly coupled, subtle in magnitude, and difficult to predict. Consequently, few computational methods have been developed and systematically benchmarked for quantitatively predicting them. Here, we establish a multiscale free-energy method and benchmark it against experimental data for a series of recently developed photoswitchable inhibitors of Escherichia coli dihydrofolate reductase (eDHFR), a crucial target in photopharmacology. Constant pH replica-exchange molecular dynamics and quantum mechanics/molecular mechanics umbrella sampling quantitatively characterize the ligands protonation-state change upon binding to the eDHFR active site. Thermodynamic integration simulations using alternative alchemical pathways, thermodynamic cycles, and protonation-state assignments were evaluated for predicting light-responsive affinity differentials and substituent effects. Direct cis-to-trans transformations with explicit treatment of environment-dependent protonation states best reproduce experimental trends. Compound-to-compound pathways are less reliable because force-field inaccuracies introduce large pK errors that are difficult to correct when protonation/deprotonation processes implicitly enter the thermodynamic cycle. TI simulations that ignore binding-induced protonation-state changes fail to consistently reproduce experimental trends. Protein-ligand and ligand-water interaction analyses further reveal the energetic and structural origins of isomer-dependent binding. This study establishes a systematic free-energy method for designing ionizable photoswitches in photopharmacology.

Antibiotic resistance in Mycobacterium tuberculosis is a pressing global health challenge demanding new therapeutic strategies. The bacterial respiratory chain comprises promising antibacterial targets, with dual inhibition of the terminal oxidases cytochrome bcc:aa3 and cytochrome bd (cyt bd) showing bactericidal activity. While bcc:aa3 inhibitors such as Q203 have advanced clinically, cyt bd remains underexplored due to difficulties in assigning activity of the purified enzyme and structurally resolving the quinol substrate binding site. Here, we report a rapid in vitro screening platform for cyt bd inhibitors by engineering a minimal respiratory system that couples the activity of cyt bd to that of a type 2 NADH dehydrogenase. This coupled assay enables spectroscopic monitoring of NADH oxidation as a proxy for cyt bd activity, allowing rapid screening of over 10,000 compounds. Screening identified WSL017, a fragment with low micromolar potency against both M. tuberculosis and E. coli cyt bd. Kinetic and structural analyses revealed competitive inhibition at the quinol-binding site, providing the first structural insights into cyt bd inhibition by a non-quinone scaffold. WSL017 displayed growth inhibition of M. tuberculosis H37ra, corroborating oxidase inhibition as a promising therapeutic strategy. This work establishes a pipeline for cyt bd inhibitor discovery and highlights new opportunities for structure-guided drug development against cytochrome bd oxidases.

Species of Bacillus bacteria including Bacillus safensis and Bacillus subtilis are finding increasing uses in biotechnology and bioremediation, thanks in part to their metabolic robustness. Malate dehydrogenase (MDH) is at the heart of central metabolism and thus a better understanding of Bacillus MDH proteins could aid in the optimization of these applications. MDH of Bacillus spp. belong to the lactate dehydrogenase (LDH)-like class of MDHs, otherwise known as the MDH3 class. Despite wide prevalence in nature among prokaryotes and archaea, this typically homotetrameric class is understudied compared to the MDH1 and MDH2 classes found in eukaryotes. We therefore recombinantly expressed and purified MDH proteins from two societally relevant Bacillus spp.-B. safensis and B. subtilis-and characterized them biophysically (via Size Exclusion Chromatography-Small Angle X-ray Scattering (SEC-SAXS) and Differential Scanning Fluorimetry (DSF)) and enzymatically (via spectroscopic activity assays). As expected based on their high sequence identity, the two MDH orthologs had similar properties in most regards, including a tetrameric structure and high susceptibility to substrate inhibition. However, we uncovered differences in conditional thermal stability, in addition to subtle differences in enzymatic activity that offer insight into the workings of LDH-like MDH. Summary statementMalate dehydrogenase (MDH) is a fundamental metabolic enzyme, from microbes to mammals, yet comparably little is known about microbial MDH, especially MDH of the tetrameric MDH3 class. We compare the biophysical and enzymatic properties of two such enzymes from the societally relevant bacterial species Bacillus subtilis and Bacillus safensis, offering useful insight with potential biotechnological implications.

Malic acid is a C4 dicarboxylic acid traditionally produced from petroleum and widely used in the food industry. As a sustainable alternative, it can also be produced as a value-added platform chemical from biomass. Previously, the Escherichia coli strain XZ658 was engineered to produce L-malate via the carbon-fixation reductive branch of the TCA cycle. In this study, we further improved this system by relieving allosteric regulation of citrate synthase, addressing redox imbalance, and enhancing malate export. These modifications approximately doubled the L-malate titer in the final strain MO128 compared to XZ658 under simple batch fermentation conditions. The process achieved a high mass yield of 1.2 g malate g-{superscript 1} glucose, highlighting the carbon-fixation capacity of the reductive TCA pathway for fermentative malate production.

Several approaches are available to increase the efficiency of CRISPR/Cas9-based genome editing, including the co-inactivation of a gene that mediates the cytotoxic activity of a compound which can be used to enrich the population in edited cells. Here we show in multiple cell lines how inactivating HSD17B11, a non-essential Short-chain Dehydrogenase/Reductase, confers a strong resistance (29- to 130-fold resistance) in both human and mouse cells to a Phenyl diAlkynylCarbinol compound (PAC) without impacting cell viability and proliferation. We show how co-inactivating HSD71B11 along with selection with PAC is usable to quickly identify efficient guide(s) against a gene of interest and to readily isolate fully inactivated clones. Altogether, this work provides an experimental framework for the facile generation of knockouts using PAC for selecting successfully inactivated cells.

Ribosomally synthesized and post-translationally modified peptides (RiPPs) are produced by biosynthetic enzymes that modify genetically encoded precursor peptide backbones and side chains. Genome mining and bioinformatics analyses targeting the multinuclear nonheme iron oxidative (MNIO) enzyme family led to the identification of a Streptomyces thermodiastaticus JCM 4840 RiPP biosynthetic gene cluster, the std cluster, which includes multiple biosynthetic enzymes and a precursor peptide containing a conserved SNKEWQE motif. Using in vitro approaches, we elucidated the modifications installed by the std biosynthetic enzymes. First, a YcaO-TfuA pair thioamidates the backbone of asparagine. Next, a peptidase S8/S53 domain fused to a NodU-like carbamoyltransferase that both carbamoylates the {varepsilon}-amino group of lysine to produce the non-proteinogenic amino acid homocitrulline and cleaves the C-terminal EWQE motif. Finally, a partner protein-MNIO pair bis-hydroxylates the {beta}- and {gamma}-carbon positions of the installed homocitrulline. The formation of homocitrulline and its subsequent modification are unprecedented in RiPP biosynthesis. Moreover, these findings expand the substrate scope of YcaO-TfuA enzymes and MNIOs and identify new roles for carbamoyl transferases in RiPP biosynthesis.

Isoflavone hydroxylases (IFHs, CYP81E) convert isoflavone aglycones into their respective hydroxylated intermediates, which direct legume isoflavones into specialized defense pathways. In soybean, their functions have been studied mostly in the context of the daidzein-derived glyceollin biosynthesis. Here we combine metabolomics-guided feature mining, phylogenetic analysis, heterologous enzymology, structural elucidation, and in planta metabolite validation to determine the functional landscape of the soybean IFH family. Analysis of a soybean isoflavonoid-enriched metabolomic dataset revealed unidentified hydroxyisoflavone features that co-accumulated with glyceollins, indicating branch chemistry that is not well-recognized. The systematic characterization of the repertoire of soybean CYP81E has demonstrated that 9 out of 11 GmIFHs are catalytically active and collectively span both 2'- and 3'- hydroxylation of the major soybean isoflavone aglycones. Among them, GmIFH9A showed broad substrate scope and regioselectivity, yielding canonical and previously unknown hydroxylated isoflavone products. NMR and LC-MS/MS were used to identify and validate the hydroxylated isoflavone products as 2'-hydroxyglycitein and 2'-hydroxyformononetin, whose presence was also confirmed in soybean roots, thus confirming two of the hidden soybean isoflavonoid network metabolites. Kinetic studies also indicated that, although the majority of GmIFHs prefer daidzein and genistein as substrates, a few isoforms are active towards methoxylated isoflavones as well, indicating functional divergence in this expanded family. Our findings collectively redefine soybean IFHs as a multi-functional enzyme module that expands the hydroxyisoflavone chemical space and reveals new biosynthetic entry points beyond canonical glyceollin pathway.